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Gerald Knizia and Garnet Kin-Lie Chan 

Department of Chemistry, Frick Laboratory, Princeton University, NJ 08544, USA 

We introduce DMET, a quantum embedding theory for computing frequency-independent quantities, such as 
ground-state properties, of infinite systems. Like dynamical mean-field theory (DMFT), DMET maps the bulk 
interacting system to a simpler impurity model and is exact in the non-interacting and atomic limits. Unlike 
DMFT, DMET is formulated in terms of the frequency-independent local density matrix, rather than the local 
Green's function. In addition, it features a finite, algebraically constructible bath of only one bath site per 
impurity site, with no bath discretization error. Frequency independence and the minimal bath make DMET 
a computationally simple and efficient method. We test the theory in the ID and 2D Hubbard models at and 
away from half-filling, and we find that compared to benchmark data, total energies, correlation functions, and 
metal-insulator transitions are well reproduced, at a tiny computational cost. 



Dynamical mean-field theory [Ij-H (DMFT) has developed 
into a powerful embedding framework for bulk quantum sys- 
tems. Its central idea is to self-consistently map the infinite 
bulk system onto an impurity model with only a few inter- 
acting impurity sites embedded in an infinite non-interacting 
bafhiQJO. In many settings, such impurity models can be 
solved using high-level many-body methods (so called impu- 
rity solvers), 16 UT41 owing to the small number of interac- 
tions. Through the bath embedding, DMFT yields predictions 
that closely approach the bulk limit, despite the greatly sim- 
plified treatment of interactions. 

The basic quantum variable in DMFT is the local Green's 
function, igjj(co) = (aJ[0J — (H — E)]~ 1 aj). As a function of 
frequency, it provides access to the local density of states as 
well as to static quantities such as energies. However, there 
are reasons to consider simpler frequency-independent quan- 
tum variables, too. For many applications frequency informa- 
tion is not required; for example, energies can be calculated 
from time-independent states alone, as can energy derivatives 
such as compressibilities or static correlation functions, and 
many other properties. Additionally, calculating stationary 
states, such as the ground-state, is much easier than calcu- 
lating Green's functions, not the least due to the practical 
availability of many powerful numerical techniques (e.g., the 
density matrix renormalization group and its tensor network 
extensions ifTTl \W\ . coupled cluster and configuration inter- 
action theories ifTTl . and lattice diffusion fljO , auxiliary field 
fl9l . and variational Monte Carlo [20 21|). An embedding 
framework based on a frequency-independent variable offers 
a potentially more efficient as well as a more flexible route to 
access static properties of bulk systems, including the possi- 
bility of using ground state methods as impurity solvers. 

Here, we propose a density matrix embedding theory 
(DMET) with the following features: (i) the infinite bulk prob- 
lem is mapped onto a self-consistent impurity problem, con- 
sisting of interacting impurity and non-interacting bath sites, 
(ii) the single-particle density matrix (ajaj) is the quantum 
variable, rather than the Green's function, and no frequency- 
dependent quantities appear in the theory, (iii) the bath repre- 
sentation consists of a single bath site per impurity site (which 
is sufficient to exactly capture embedding effects at the mean- 
field level) and (iv) the bath can be constructed algebraically 



without any fitting. Feature (i) is analogous to DMFT, and 
as we show below, the basic physics of DMET is similar to 
DMFT. Features (ii)-(iv), however, are different. They allow 
for the primary numerical advantage of DMET: computing 
ground-state properties of a cluster model with L impurity 
sites requires only solving for the ground-state of a cluster 
plus bath problem of size 2L, and this is much cheaper than 
the corresponding DMFT calculation. 

To motivate the DMET construction, we first consider an 
exact single-site embedding of the infinite lattice Hubbard 
model, with Hamiltonian 

H= £ ta] a a j(y + Y,U n^nn . ( 1 ) 

For simplicity, we focus on ground-state properties, but we 
outline the extension to excited and thermal states below. The 
ground-state I 1 ?) of H can be mapped onto the ground state of 
a simple impurity model, consisting of a single impurity site 
embedded with a single bath site and with Hamiltonian H'. 
This follows from the Schmidt decomposition of I 1 ?), I 1 ?) = 
Y,i A; 1 a,-) |j3,-) l22l . where |a,) are states of a single Hubbard 
site (viewed as an "impurity") and |/3,) are states in the Hilbert 
space of the remaining lattice sites. Note that the number of 
1/3,) states, M, equals the number of states of the single impu- 
rity site. Consequently, |j3,-) can be interpreted as the states 
of a single bath site. With this identification, the exact impu- 
rity Hamiltonian H' (with the same ground-state as H) may be 
constructed by projecting H onto the Schmidt basis of impu- 
rity and bath states, H' = l^iPj) ( a i'A/|#|a;'j3/) («/'j3/|- 
More explicitly, H' can be expressed using the impurity and 
bath fermionic operators av\ 

H' = c + £ epblbo + v(a T a b a + b T a a a ) 
a 

+ v'm a (a a ,b a , +b a ,a a ,) + Un t n l + Upm t m l (2) 

where n a — a' a a a , m a = b]jb a , and c, Ep , v, v' , Up follow from 
matrix elements of H with the Schmidt basis, for example, 
V = (aoj3o|//|o;ij3i). Because the impurity Hamiltonian H' 
has the same ground-state I*?) as the Hubbard Hamiltonian 
H, expectation values of the infinite lattice can be obtained 
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exactly from the impurity model. (This is very different from 
DMFT, where expectation values of the lattice and impurity 
model are not related in a simple way). For example, the local 
density matrix is given by (a~ a a a ) H i = (fig ao CT )# (the latter 
referring to lattice sites, eq. ([T|), while the lattice Hubbard 
energy per site, E, is obtained from the terms in H' which 
contain the impurity ay> operator, E = v(a' a b a + b a a a ) + 

La^a' v ' m <y {a a >b& + b a ,a & ) + U n t "|)//' ■ 

The above shows that an impurity model with a single im- 
purity site and a single bath site is in principle sufficient to 
exactly represent ground state properites. However, the ex- 
act construction is not practically useful: the bath terms in H' 
require knowledge of the interacting ground-state I^P) on the 
infinite lattice and its Schmidt decomposition. The basic idea 
in DMET is to replace the exact embedding of the Hamilto- 
nian H by one that is exact for a one-particle mean-field lat- 
tice Hamiltonian h. The corresponding mean-field embedding 
bath terms are then easy to compute because the ground-state 
of h is a Slater determinant |<t>), and its Schmidt decompo- 
sition is easily obtained, at a cost no greater than the one- 
particle diagonalization of h itself [22—24]. As h we choose 
the one-particle part of H combined with an on-site mean-field 
interaction operator u (to be determined), 

h = E ta \o a ie +Y, u ( n n + n il)- ( 3 ) 
(ij)o io 

The corresponding DMET impurity Hamiltonian i/i mp con- 
tains a single interacting impurity site, now embedded with a 
non-interacting bath site, 

#imp = E (v{alb a + bla a ) + ub r a b a ) + C/n t n;. (4) 

a 

The terms involving the bath operators b^> in flj mp are con- 
structed analogously to the bath terms in eq. Q, by pro- 
jecting h onto the Schmidt basis of |<t>), e.g. given |<J>) = 

Ylf Af| a,-) |J3,-}, v = (aoj3i|/i|aiA)). Note tnat the batn terms 
arising in this way contain only one-particle operators. The 
mean-field interaction operator u, which defines both the lat- 
tice Hamiltonian h and (indirectly) the impurity Hamiltonian 
i/imp, is analogous to the DMFT self-energy. By changing u, 
we also change the bath terms in H\ m? , and we can adjust u 
to obtain a mean-field embedding that "optimally" mimics the 
exact embedding. One way is to require full self-consistency 
of the local density matrix, that is, that the mean-field lat- 
tice Hamiltonian h and the DMET impurity Hamiltonian H imp 

yield the same local density matrix: (affa a )H imp = {a\ a a.Qa)h- 
This is similar to the self-consistency condition on the Green's 
function in DMFT. We here use a slightly different self- 
consistency condition, which we have found to be numerically 
favorable. We minimize the difference between the total den- 
sity matrices (using both impurity and bath operators) evalu- 
ated for the ground states of h and //; mp , 

mm £ ((4 c cT>H iffl p - <4 c <r>ft) • (5) 
o,ce{a,b} 



This approximately maximizes the overlap between the mean 
field wave function and the full wave function in the impu- 
rity model. Other choices are possible, similar to the differ- 
ent choices of self-consistency condition in DMFT and self- 
energy functional theories J5] [25) . 

In the exact embedding construction, the impurity Hamil- 
tonian H' and the infinite Hubbard Hamiltonian H shared the 
same ground-state, thus exact expectation values could be ob- 
tained from the impurity model. In the case of the DMET im- 
purity Hamiltonian //i mp , this is no longer the case, because 
the embedding is constructed for the mean-field h rather than 
H. Nonetheless, in DMET we make the approximation that 
lattice quantities are approximated by the expectation values 
of the impurity model. The local density matrix in DMET is 
thus defined as {a^ a a a )H imp , and the energy-per-site is 

E =Y, v ( a c t, G+ b ° a a)H nw +U(n t n l ) Himp . (6) 

a 

So far we have discussed single-site DMET. Cluster exten- 
sions can also be formulated. The simplest cluster extension 
of DMET for L impurity sites is analogous to L site cellular 
DMFT 0, and the impurity Hamiltonian //; mp is given by 

tfimp = E E \ v iM]a b io +b] a a ja ) +Uijb] a b ja \ +Un nnil 
a i,j=l 

(7) 

where the bath , interaction uu, and coupling Vu opera- 
tors are all generalized to L sites. The construction of the 
couplings from the Schmidt basis of h, the self-consistency 
conditions on u, and the approximation of lattice quantities by 
impurity quantities, are all directly analogous to the single-site 
case. Note that this kind of cluster DMET breaks translational 
invariance, similarly to cellular DMFT [5 1. 

Further extensions are possible, but not followed here. 
While we have focused on ground-states, the DMET construc- 
tion is equally applicable to excited states (by carrying out the 
Schmidt decomposition for the excited state, rather than the 
ground-state). In the case of thermal states, we would regard 
exp(— £>H) as a state vector in the enlarged Liouville space to 
carry out the Schmidt decomposition, leading to impurity and 
bath sites with twice the number of degrees as in the case of 
pure states. 

How does DMET perform? Similarly to DMFT, the DMET 
construction is exact in the non-interacting (U = 0) and atomic 
(t = 0) limits of the Hubbard model [3 ]. In the non-interacting 
limit, the mean-field lattice Hamiltonian and Hubbard Hamil- 
tonian are the same, h = H, thus the mean-field embedding 
used in DMET is exact. In the atomic limit, the Hubbard 
Hamiltonian decouples the sites, the impurity-bath coupling 
v = vanishes, and the impurity expectation values are iden- 
tical to Hubbard model expectation values. DMET thus acts 
similarly to DMFT in providing an interpolation between 
metallic and Mott insulating behaviour. Note also that DMET 
has the same diagrammatic structure as DMFT; solution of 
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FIG. 1. (color online) E(n) with various cluster sizes: (a) ID Hubbard model, and (b) 2D Hubbard model. In both cases good agreement with 
reference data is found. For (b), QMC reference data is taken from Refs.l26landl27l 
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FIG. 2. (color online) ID Hubbard model: Local double occupancy 
(n 1 ^ ■ n]) as function of U; at half-, 2/6, and 1 /6 filling. Except for 
the 1-site DMET case, the reference shapes are well reproduced. 



FIG. 3. (color online) Metal-insulator transition in the lD-Hubbard 
model: a Mott gap is seen for cluster sizes > 1. The size of the gap 
is accurately reproduced. 



the impurity Hamiltonian with an exact local interaction en- 
sures that the resulting observables have diagrams where the 
local interactions are treated to all orders. The different self- 
consistency condition and bath construction, however, mean 
that the theories are not identical even in the static limit. 

To assess the validity and accuracy of DMET away from ex- 
act limits, we now turn to numerical DMET studies of the ID 
and 2D Hubbard models, as a function of U and filling. These 
models were chosen because high-quality reference data is 
available. In ID, we compare to exact results from the Bethe 
ansatz[28, 29], while in 2D we compare to recent auxiliary 
field quantum Monte Carlo (QMC) calculations ESl |27l . In 
the numerical DMET calculations, the infinite lattice used to 
define the mean-field Hamiltonian h in eq. ([3]) is approximated 
by a large lattice with antiperiodic boundary conditions, us- 
ing 480 sites in ID case, and 24x48 sites in the 2D case. 
The resulting finite size errors are substantially smaller than 
the intrinsic errors from the DMET approximation. For the 
DMET impurity solver we employed an exact diagonalization 
algorithm OUl . The full source code for the prototype DMET 
and the impurity solver is freely available ll24l . 

We first discuss DMET's ability to describe energetics and 
correlation functions in the ID and 2D Hubbard models. 
Fig. [T] gives the energy per site E, as a function of site occu- 
pancy n, using DMET of various cluster sizes. Note that while 



the ID Bethe ansatz results are exact, the 2D AFQMC refer- 
ence data is not [26|, except at half-filling ll27l . As seen from 
the figure, the overall shapes of the reference E[n) curves are 
reproduced well. Unsurprisingly, DMET is most accurate at 
smaller U. The accuracy increases with impurity cluster size, 
but even single-site DMET produces reasonable energies: at 
half-filling in ID, single-site DMET is accurate to 0.12% at 
U = 1, and 4% at U = 4, while at half-filling in 2D, it is ac- 
curate to 1.8% at U = 4. In Fig. [2] we plot the local double- 
occupancy («|«|), which is a measure of the Mott insulating 
character of the state [31|. DMET also gives a good descrip- 
tion of this quantity, although at half-filling, single-site DMET 
does not capture the correct change in curvature of as 
a function of U due to the neglect of short-range singlet for- 
mation, while larger clusters recover this interaction effect. 

A central question in the DMFT treatment of the ID and 
2D Hubbard models is the location and nature of the metal- 
insulator transition and the corresponding size of the Mott 
gap Il32l - l34ll . In the antiferromagnetic DMFT solution a gap 
opens at U = due to perfect nesting. The paramagnetic 
DMFT solution is more complicated, and there is a critical 
interaction strength U, which depends on the cluster size, 
beyond which a Mott insulator paramagnetic solution is sta- 
bilised. In a ground-state calculation, the opening of a Mott 
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FIG. 4. (color online) Comparison of DMET and DMFT descriptions 
of the Metal-insulator transition. DMFT data is taken from Ref.l36l 



gap can be detected by a vanishing compressibility, d(n) /djj,. 
In Fig. [3]we plot (n) against ji for DMET calculations on the 
ID Hubbard model. In the antiferromagnetic DMET calcu- 
lations, a gap opens at U = as in DMFT. In the paramag- 
netic DMET calculations, a Mott transition is observed in the 
cluster DMET calculations for larger U's, although not in the 
single-site DMET. Note that single-site DMFT also does not 
display a Mott transition at moderate U l35l . Compared to 
the Bethe ansatz, the Mott gap is well reproduced, with the 
error decreasing with cluster size. The n(ji) behavior close 
to half-filling shows a bi-stability, indicating that the metal- 
insulator transition is likely to be first-order, similar to what 
is found in DMFT. As seen from Fig. |4] the shape of the 
n(jj.) curves resembles that obtained in recent cellular DMFT 
studies 1 35 36 1. In the 2D Hubbard model, using a 2 x 2 clus- 
ter DMET, we find co-existence of metallic and insulating 
paramagnetic solutions over a range of U, in addition to the 
antiferromagnetic solution (which is correctly identified as the 
lowest energy phase). The co-existence region starts at about 
U « 6 and extends to U ~ 10.5. However, the metallic solu- 
tion is clearly favored up to U ss 9.5 as seen from Fig. |5](left). 
Cluster DMFT [32, 33 1 and variational plaquette l34l calcula- 
tions also show a significant co-existence region for paramag- 
netic metallic and insulating solutions with the metal-insulator 
transition occurring near U ~ 6. In Fig. |5](right) we plot n(jx) 
for U = 12. We see that the structure away from half-filling is 
more complicated than in the ID plots. Note that this behavior 
occurs at moderate doping where complex phases of the Hub- 
bard model are expected to exist. Overall, we find that the 
DMET and DMFT descriptions of the paramagnetic metal- 
insulator transition are similar, but the DMET calculations of 
ground-state properties are much cheaper: each one requires 
computing only the ground-state on a very small number of 
sites, which takes only seconds. 

To summarize, here we introduced a density matrix em- 
bedding theory (DMET) to obtain properties of infinite bulk 
interacting systems from a simpler embedded quantum impu- 
rity model. While similar to DMFT, DMET is formulated in 
terms of the frequency independent local density matrix. The 
absence of frequency yields both formal and practical advan- 



FIG. 5. (color online) 2D Hubbard model: (left) Phase stability 
of metallic and insulating solutions at half filling as function of U ; 
(right) Paramagnetic metal-insulator transition at U = 12. 



tages. Further, the DMET bath construction requires only a 
small number of sites (one per impurity), and is obtained al- 
gebraically without non-linear fitting. There is no bath dis- 
cretization error. The simple bath construction and lack of 
frequency means that DMET calculations for static properties 
are much simpler and faster than the corresponding DMFT 
calculations. We showed that DMET contains similar local 
physics and is exact in the same limits as DMFT, namely for 
weak interactions and weak couplings. Furthermore, for the 
ID and 2D Hubbard models, accurate behaviour for the en- 
ergetics, correlations, and metal-insulator transition was ob- 
served. Consequently, we conclude that DMET provides an 
appealing alternative to DMFT for the ground-state properties 
of infinite interacting systems. Applications to further prob- 
lems, as well as extensions e.g. to long-range interactions, are 
now underway. 
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